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We review the status of non-perturbative analyses of multi-jet event shape distributions and mean values, highlighting 
the physical insight on QCD dynamics they can provide. 

1. INTRODUCTION 

In every QCD observable, perturbative (PT) and non-perturbative (NP) dynamics are inseparable, the reason being 
that the observed degrees of freedom are not quarks and gluons, the elementary particles entering the QCD bare 
Lagrangian, but hadrons, whose description in terms of partons is well beyond the domain of PT theory. Fortunately, 
at least for sufficiently inclusive observables, the difference between parton and hadron level predictions is suppressed 
by inverse powers of the hard scale of the process. As far as the latter is short-distance dominated, one can safely 
compute QCD observables using the PT parton language, and interpret the discrepancy with experimental data, in 
case any is seen, as the need for NP hadronisation corrections. One can then adopt two complementary approaches. 
One is to consider observables whose hadronisation corrections are almost negligible, for instance total cross sections 
or inclusive non-QCD particle distributions. These observables can be computed in PT QCD and exploited to 
determine the value of the coupling a s 0] . The other approach is to consider observables which are very sensitive to 
NP physics, in order to have an insight on the hadronisation mechanism. The best known example is event shape 
variable distributions and mean values. These variables are constructed by combining final state momenta to obtain 
a number that gives an idea of the geometrical properties of hadron energy-momentum flow. The value of an event 
shape V is related to the scale at which hadrons are probed, so that measuring event shape distributions makes it 
possible to study physics at very different scales, which range from the domain of PT QCD (V ~ 1) down to the 
confinement region (V ~ Aqcd) where the quark/gluon language is scarcely applicable. For this reason, although 
shape variables have been used to measure a s , they are ideal for investigating properties of QCD dynamics. So 
far, both experimental and theoretical investigations have been restricted only to two-jet event shapes (see |3] for a 
review). Here we will discuss how the extension of such studies to multi-jet event shapes can shed further light on 
the interplay between PT and NP physics in QCD observables. 



2. TWO- JET STUDIES 

2.1. NP correction to event shapes 

Experimental data clearly indicates that PT QCD alone is not enough to predict event-shape distributions and 
mean values Let us consider the mother of all event shapes, the thrust 0] 

T = max E ^;" T| , r = l-T. (1) 

T,h \Ph\ 

The thrust is a measure of particle alignment, and is a typical two-jet variable, since it vanishes in the limit of two 
narrow jets. As one can see from fig.^ the dependence of the mean value of 1 — T on the e + e _ centre-of-mass energy 
Q is correctly described only after adding to the QCD fixed order prediction a NP 1 / Q-suppressed correction 



>pt + ( t )np' ( t )pt = a s(Q) T i +& S (Q) t 2 + . . . , (r) NP ~ — , (2) 



Cr 

Q 



where C T ~ lGeV when (r) PT is evaluated at next-to-leading order (NLO). 
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Figure 1: The mean value of 1— T as a function of Q = y/a. Experimental data from various collaborations are plotted against 
theoretical predictions. The figure is taken from ref. [jj . 



One might think that such a discrepancy could be removed by including higher orders in the PT expansion. Actu- 
ally Sterman observed that a term 18 can already mimic al/Q behaviour |5|. However, theoretical analyses show 
that the PT series is an intrinsically ill-defined object, since it is doomed to diverge factorially (for a review, see Q). 
Attempts to regularise such a divergence give rise to a power-suppressed ambiguity, known as infrared (IR) renor- 
malon. Looking more closely at the origin of the divergence, one can see that it arises when resumming renormalon 
chain graphs containing an arbitrary number of linked quark or gluon bubbles Q]. A renormalon resummation for 
the thrust mean value yields 

AC F a s ^ ( 2(3 a s \ n f dk t a s (k t ) k t i , 

(r) PT * — - gn! j x 2C F J —drj — -e , (3) 

where /3q and fti are the first two coefficients of the QCD beta function. After a renormalon analysis, the series in 
eq. J2J gives a 1/Q ambiguity. The last equality in eq. © results from the fact that the series can be seen as the 
PT expansion of the integral of the running coupling a s (k t ) down to the infrared. Since the PT coupling develops 
a Landau singularity at low momenta, one may be tempted to ascribe the divergence to the presence of the Landau 
pole in the k t integration contour. However, eq. shows that the divergence of (t) pt is determined only by /3rj and 
/3i, i.e. it is independent of the particular coupling adopted. This naive observation is confirmed by more refined 
theoretical analyses, which show that IR renormalons are always present whatever is the behaviour of the coupling 
at low scales Q. 

The main message of this discussion is that, in order to obtain a satisfactory theoretical understanding of event 
shape observables, PT theory must always be supplemented with information on QCD dynamics at low scales. In 
particular, a size of the power correction of around lGeV suggests that NP effects arise from partons which have 
started the blanching process that leads to the formation of hadrons. 



2.2. Two-jet shapes and the Feynman tube model 

Consider now the specific case of an event-shape V in e + e~ annihilation that vanishes in the two-jet limit (thrust, 
C-parameter, heavy-jet mass, etc.). The mean value (V) will receive its main contribution from hard particles whose 
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transverse momenta (with respect to the thrust axis) are of the order of the hard scale Q. There are however soft 
hadrons, with transverse momenta up to about lGeV, whose contribution to (V) cannot be fully computed with PT 
techniques. The Feynman tube model Q gives a phenomenological description of such hadrons by assuming that 
their distribution is uniform in rapidity: 

dn h 



din ktdn 

Observing also that their contribution 6V to V is additive 



■ (4) 



Q 

one obtains a l/Q-suppressed NP correction (V) NP to (V), given by 



W-E^M*). (5) 



h 

where cy is a calculable coefficient and 



WNP=E/^^(^)|M^) = ^^, r y =\ dij J\ (//! . ((,) 



( fc *>NP = E / ^h(h)h, (7) 



is the average transverse momentum of the produced hadrons, and represents a genuine NP quantity. 

The assumption of eq. Q makes it possible to completely factorise rapidity and transverse momentum dependence. 
Furthermore, eq. © implies that power corrections to two-jet event shapes are universal, in the sense that, after 
computing the variable dependent coefficient cy, they depend only on the NP parameter (fct) NP , which is the same 
for all variables. 

A comment is in order concerning the validity of eq. Since in a hard collision most particles are aligned along 
the thrust axis, one expects that the distribution of hadrons away from the jets should be invariant over small boosts 
along the thrust axis direction, but that, as soon as one moves forward in rapidity, eq. gets modified. Fortunately, 
for most event shapes fv(v) strongly suppresses the contribution of hadrons at large rapidities, so that leading power 
corrections are determined only by soft hadrons in a central rapidity region, which are well described by eq. (Q. 1 

This very same approach can be used to compute power correction to event-shape distributions. Since V — Vpt+5V 
with SV ~ e/Q, the distribution da/dV can be seen as a convolution of a PT distribution dapx/dVpT and a NP 
shape function /np(e) PJJ 

~% = I dVpT de -^(^t) / Np ( £ ) 6 ( V - VpT - ?0 • ( 8 ) 
o- dV J a dVpT \ Q J 

In the region VQ 3> (e) one can expand /np(c) around (e) and obtain 

1 da 1 da PT fy_(±\ (Q) 



adV a dV \ Q 

Recalling now that (e)/Q — (V) NP , one observes that leading power corrections to da/dV result simply in a shift of 
the corresponding PT distribution, whose magnitude is the same as the power correction to (V) 

To conclude, the basic assumption of the Feynman tube model, that the distribution of soft (central) hadrons is 
uniform in rapidity, implies that leading 1/Q NP corrections both to mean values and distributions can be expressed 
in terms of a single universal parameter (fct) NP . This universality hypothesis has been tested against experimental 
data and is found to hold within 20% accuracy 0. 



1 Surprisingly enough, eq. J1J gives the correct leading power correction also for variables that are uniform in rapidity. This is due to 
a tricky interplay between PT and NP QCD radiation, which will be discussed in the next section. 
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3. EXTENSION TO MULTI-JET EVENT SHAPES 
3.1. Intra-jet hadron distribution 

When moving from two-jet to multi-jet events one encounters difficulties in extending eq. ijlf. One can reasonably 
think that particles inside one of the jets (inter-jet hadrons) are produced uniformly in rapidity and azimuth with 
respect to the jet axis. However, for particles at large angles with respect to all jets (intra-jet hadrons), since there 
is no natural way to define r\ and <fi, one does not expect dn^ to have a simple expression. To find a reasonable 
assumption for dn^, we start from two-jet events and consider the PT probability dw(k) for the emission of a soft 
dressed gluon k off a back-to-back quark- antiquark system (whose momenta are p and p) in a colour singlet 

a OA n dk t a d^OsiMt) 1, pk 2 (2pk)(2kp) 

dw (k)=C FTf dn-—, r,= -ln-, k t = — , (10) 

where n and k t , the gluon rapidity and transverse momentum, have been written in a Lorentz invariant form and the 
coupling is taken in the physical CMW scheme 01 . The inter-jet hadron distribution in eq. can be obtained from 
eq. Hl()|) by replacing 2 CFa s (k t )/Tr with &h(k t ). In this way dny t might be interpreted as an effective measure 
of QCD interaction strength in the infrared. When the number of hard emitters is larger than two, the distribution 
dw(k) can be written as a sum of contributions from all dipoles formed by the hard partons [l^j : 

where now Kij, rjij and (/>ij are the transverse momentum, rapidity and azimuth in the (ij)-dipole centre-of-mass 
frame, and the coupling is again in the CMW scheme. The vector Tj represents the colour charge of hard parton 
Pi. From colour conservation one has that for less than four emitters Ti ■ Tj are numbers, while starting from four 
partons they are actual matrices in colour space. 

We now postulate that, after a suitable extension of the CMW coupling at low scales has been introduced, eq. 
gives the distribution of intra-jet hadrons in a multi-jet event. This corresponds to the so-called local parton-hadron 
duality (LPHD) hypothesis (see ^| for a recent review), which states that hadron flow is determined by parton 
flow. This assumption has solid phenomenological bases, since it describes very well hadron multiplicity flows |lq 
and string/drag effects in three-jet events [r?j |. 

Before constructing such an extended coupling, we first discuss to what extent the soft large-angle approximation 
is sufficient to compute 1/Q power corrections to event shapes. In fact, Kij can reach the NP domain not only for 
soft large-angle, but also for hard collinear emissions. Considering emissions collinear to a given leg, we can classify 
event shapes according to whether they are damped in rapidity or not. 

If V is damped in rapidity, soft collinear particles with transverse momentum ku and rapidity rji with respect to 
the given leg contribute to V with a correction 

^-E^MW), fvW-e-", V^oo. (12) 
i y 

This implies that only large angle emissions contribute to SV, up to higher power corrections. 

If V is uniform in rapidity, we have SV ~ J^. hti, so that one might think that all rapidities contribute equally to 
SV. However, due to PT radiation, each hard parton takes a recoil, and acquires a transverse momentum p t ~ VQ. 
But one measures transverse momenta with respect to a fixed axis (e.g. the thrust axis), and not with respect to 
the hard emitter's direction (the jet broadenings Q are a relevant example). Since, as we have seen in the previous 
section, leading power corrections come from the region SV -C V, we must have then ku/Q ~ SV -C Pt/Q ~ V. 
This means that particles giving the leading contribution to SV must be soft, and displaced from the hard emitting 
parton, i.e. soft at large angles. To conclude, eq. i|ll[l gives the physically correct starting point to compute 1/Q 
power corrections to event-shape variables. 
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3.2. Dokshitzer-Marchesini-Webber (DMW) extension of CMW coupling 

The CMW coupling can be extended at low scales via the dispersion relation 

«-(«) = - / jJr^) ^™ 2 ) = « 2 I { J™ k2) 2 <W™ 2 ) > (13) 

where a e g(m 2 ) is the logarithmic derivative of the spectral density p s (m 2 ), and is itself a QCD coupling. 

This approach automatically incorporates the LPHD hypothesis by assuming that hadronisation corrections are 
due to gluers, extra-soft gluons with transverse momenta k ~ Aqcd , whose emission probability is ruled by the NP 
part of the dispersive coupling. 

Within the dispersive DMW approach, power corrections to event shapes can be computed by considering the 
emission of a single gluer, and following a standard procedure [2p|: 

1. real and virtual corrections are combined to obtain the dispersive representation of the QCD coupling. This 
means that in eq. (|llf> . for each dipole, one can make the substitution 

"»(«) f°° dm2 , 2n /,„>> 

r a c ft(m ) ; (14) 



Iq (to 2 + K 2 ) 2 

2. the gluer is given a mass to 2 , i.e. it is allowed to decay inclusively. In practice this results in replacing 
everywhere in the event-shape definition k 2 with k 2 + to 2 , so that the contribution to a variable SV^l T of a 
gluer emitted off the (ij)-dipole is 



<£ = ^^#W), d5) 



ill' 



where k, n and 4> are, as usual, considered in the emitting dipole centre-of-mass frame. Then the leading power 
correction to (V) becomes 



{V) »f° = . f,) , = I dn^-fPin, 0) , (16) 

where 



dn 2 dm 2 /-r, a c s(m 2 ) 2 f°° dm 2 



\/ k 2 + to 2 ~— / — —to Sa c ff(m 2 ) , (17) 



(m 2 + K 2 ) 2 V 7T 7T 

turns out to be related to the 1/2-moment of (5a e ff(TO 2 ), the NP part of the dispersive coupling; 

3. non-inclusivcncss of event shapes with respect to secondary gluon decay is accounted for by replacing (ft) naive 
with («) NP = M («) na i ve5 where M is the so-called Milan factor [lj,|21j. It is also customary to rewrite (k) np 
in terms of ao (///), the average of the dispersive coupling below the merging scale /ij, in such a way that the 
sum of PT and NP contributions to (V) is free of IR renormalons: 

4ur f^ 1 dk 

(K) NF = J i-M(ao((Jii)-a s (Q) + 0(a 2 )) , ao(pi) = —a s (k). (18) 

We are now able to discuss what kind of information on QCD dynamics we obtain from power corrections to two- 
or multi-jet event shapes. 

From eq. 1)16(1 the power correction to a two-jet variable can be written in the form 

/t/\ _ n ( k )np „ _ (^)np „ f 1Q \ 
(1/) NP - Of c v — — g — cy, (19) 

where (fct) NP is the NP parameter appearing in eq. J7J, and the last equality has been obtained by comparing eq. JBJ) 
and eq. H16[). We see then that both the dispersive approach and the Feynman tube model give the same result for 
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(V) NP . In fact universality of leading NP corrections to two-jet event shapes implies simply that central soft hadrons 
are produced uniformly in rapidity. 

When increasing the number of jets, the NP correction to a multi-jet variable is given by 

(V) NP = { -^c v , = £(-f, • fj) . (20) 

The main features of this PT QCD inspired correction are: 

1. the same parameter (k) np determines the magnitude of power corrections for all event shapes and for any 
number of hard jets; 

2. (V) NP depends on the colour charges of emitting partons through the factors T ■ Tj\ 

3. the power correction is sensitive to the geometry of the underlying hard event (the angles between the jets) 

„(*J 
-v 



through the coefficients c. 



This colour and geometry dependence is a highly non-trivial property that, once established, would imply that 
hadronisation preserves the main characteristics of partonic energy-momentum flow. 

3.3. Three-jet event shapes in e + e~ annihilation 

We now consider the specific case of three-jet events in e + e~~ annihilation. At Born level a three-jet event is made 
up of a quark p q , an antiquark pq and a gluon p g . Transverse momenta are, as usual, defined with respect to the 
thrust axis Ht, and one defines an event-plane as the one containing Ht and thrust major axis um, defined as the 
one maximising the projection of transverse momenta. Every transverse momentum k t can be decomposed into an 
in-plane component fc m and an out-of-plane component fc out as follows 

k t = fc m n hl + k out ri out , n in = ri M , "out = n T x ri M • (21) 



The two variables for which there exists a prediction for leading power corrections are the /^-parameter |22j and the 
thrust-minor T m (a.k.a. -Kout) l^ j- 

The D-parameter is defined as the determinant of the tensor 

a 

6> a/3 Q = y^f, D = 27det(9 = 27AiA 2 A 3 . (22) 



The soft particle contribution SD is given by 



K7 sin 



SD ~ 27AiA 2 > 1 yt , (23) 

where is the transverse momentum in the emitting dipole centre-of-mass frame and <pi = for an emission inside 
the event plane. Due to the presence of the energy u>i in each denominator of eq. I|23|l . this variable is damped in 
rapidity, and the coefficient cd is given by 

cd = C F c q D (T, T M ) + C F cl(T, T M ) + C A c 9 D (T, T M ) , (24) 

where c q D , c q D and c^, depend on the event geometry through T and Tm- 
The thrust-minor T m measures the out-of-event-plane momentum flow: 

h i 
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This variable is uniform in rapidity, so that one needs to consider the fact that each hard emitting parton recoils 
against PT radiation. This results in 

e Jf „, M | sto , l) ( &h ^ + Crh ^ T + c , b ^|| T ) , fl-MD-fgl-M-!, (26) 

where Q 2 ^ = 2pj ■ pj . We notice that for each hard emitter the power correction is proportional to its colour charge 
and to the rapidity interval available to the emitted gluon. The fact that the NP correction depends on hard parton 
recoil implies that the final answer will be obtained after averaging cx out over the PT recoil distribution, which results 
in a complicated but very interesting interplay between PT and NP physics. 

Fig. |21 shows the comparison between theoretical predictions [2^] and experimental data for the /^-parameter 
differential distribution and mean values. Three-jet events are selected by imposing a lower limit y cut = 0.1 on the 
(Durham) three-jet resolution j/3. The distribution in fig. [21 represents the first theoretical prediction for a three-jet 
variable at the state-of-the-art accuracy, which includes all-order PT resummation of logarithmic enhanced contribu- 
tion to next-to-leading logarithmic (NLL) accuracy, matching to the corresponding NLO distribution obtained with 
NLOJET++ [2^, and l/Q NP power corrections. For each Born configuration, the leading NP correction results in a 
shift of the PT distribution. The presented distribution is also slightly squeezed since the shift depends on the event 
geometry and one has to integrate over different Born configurations. The comparison to data is overall quite good. 
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Figure 2: Theoretical predictions [22j for distribution (left) and mean values (right) of the D-parameter compared to ALEPH 
data [3, for Q = 91.2GeV. Mean values for the data are computed from the corresponding distributions, with a 2% estimated 
experimental error and an extra 1% systematic error to account for the extraction of mean values from distributions [2t|. 



One can observe that the shift is quite large, about twice as much as in two-jet event shapes. This is mainly due 
to NP radiation off a gluon, which has a large colour charge. Although a full analysis of theoretical uncertainties 
remains to be performed [27| , the magnitude of the shift suggests that leading power corrections might not describe 
with sufficient accuracy the distribution around the peak, but that higher powers, or even a shape function, might 
be needed. 



3.4. Three-jet event shapes in hadron-hadron collisions 

One can construct three-jet event shapes also in hadron-hadron collisions. The simplest case is to consider a 
process with an electro-weak vector boson Z° or recoiling against a high-p t jet. The process involves three hard 
emitters, two incoming (p\ and P2) and one outgoing (733). 
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The event plane is the one containing the beam and the momentum q of the vector boson, so that we have 
"out = (q x Pi)/(l < zl|Pi|)- Since one does not measure inside the beam pipe, any measured hadron p^ will have 
\Vh\ < ??o, where 770 is the maximum available rapidity. Therefore the out-of-plane momentum flow -Kout has to be 
defined as 

K out = J2\Ph ut \ xefro-M), (27) 

h 

In the soft limit <5-Ko U t gets two contributions, one from particles emitted directly from the three hard partons, the 
other from hadrons produced in beam remnant interactions, the so-called soft underlying event [2^ : 

SK out = SKfr + SK—* , SKfr * E' K * I sin & I > SK~ nt * E' fe « I sin & I > ( 28 ) 

i i 

where each primed sum is restricted to particles outside the beam pipe. In the term 5K^ eT the transverse momentum 
Ki is considered as usual in the emitting dipole centre-of-mass frame, while in 5K^ nant each ku is in the laboratory 
frame. 

As before, the gluer contribution depends on the same parameter (k) np as for e + e~ event shapes, and is given by: 

5K^ = («> N p c Kmt , c KmA c (| sin 0|) (c^o - ry 3 ) + C 2 ( m + m ) + C 3 In , (29) 

where Qt is the transverse momentum of the vector boson, which represents also the hard scale of the process, and 
Ci is the colour charge of hard parton pi. As in the e + e~ case, for each emitting hard parton, the power correction 
is proportional to the rapidity interval available to the emitted gluer. 

The beam remnant contribution can be computed by assuming that the distribution of hadrons produced in soft 
beam-remnant interactions is uniform in rapidity, which gives 

8K~ nt = (I sin 0|) x (^) rcmnant x 2r l0 . (30) 

We have then access to the new NP parameter (fc t } remllarlt , the average transverse momentum of hadrons produced in 
beam remnant interactions. This very same parameter will appear in power corrections to event shapes in hadronic 



dijet production [3j. 



4. A RICH PHENOMENOLOGICAL PROGRAM 

The considerations we have presented so far can be seen as the starting point for a rich phenomcnological program. 
First of all one has to keep in mind that multi-jet event-shape observables are not ideal to perform pure PT studies 
aimed at precise determinations of the strong coupling a s . Large theoretical uncertainties arc associated to PT 
higher orders and to NP power corrections, and jet selection cuts reduce the number of available events. However, 
from the point of view of understanding QCD dynamics, they are an incredibly valuable tool. A research program 
has already started with the study of power corrections to three-jet event shapes in e + e~ annihilation and DIS |27|. 
aimed at determining whether the NP parameter «o is the same as for two-jet event shapes. Once the universality of 
ao is established, one could move to considering hadron-hadron collisions. There one could test our understanding 
of beam-remnant interactions in a complementary way with respect to what has been done so far. At the moment 
the only observable used to test models of the underlying event is the away-from-jet particle/energy flow |3o| |. for 
which we have a less clear theoretical understanding than we have for event shapes. Provided our modelling of the 
underlying event proves to be correct, we have access to the NP parameter (fct) romnant , and are able to use this 
parameter to study event-shapes in hadronic dijet production |3l|. where a lot of progress has been done in recent 
years, especially since when the automated resummation program CAESAR has become available 
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